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ABSTRACT 


A theoretical model for the prediction of the temperature 
field and the concentration field of diffusing species in the 
vicinity of the cutting edge during orthogonal machining has been 
developed. The finite element method has been applied to solve 
the steady heat transfer and transient species transport equations 
in two dimensions. The changes in the thermal properties and 
diffusion coefficient with temperature have been accounted for. 
Based on the temperature and concentration solutions, changes in 
the strength of the tool material with time at each location have 
been predicted, assuming carbon diffusion. Comparing the local 
strength with the actual stress, regions which are likely to 
undergo wear have been identified. A few typical machining 
situations have been analysed and results show reasonable 
agreement where comparisons are possible. The present approach 
has promise for being a valuable theoretical aid for analyzing 
tool wear due to diffusion at high cutting speeds. 
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CHAPTER 1 


INTRODUCTION 


1.1 GENERAL BACKGROUND 


The importance of machining in the modern engineering 
industry need not be over-emphasised when it is stated that nearly 
half of the engineering products are produced through machining at 
one stage or the other. In the age of wide spread automation, we 
are witnessing the trend that the demand for productivity and 
process optimi2ation is increasing day by day. As the search for 
new materials and processes is continuing, effort is still 
underway to improve conventional process by reducing the material 
wcistage to a bare minimum, increasing tool-life and improving the 
product finish. In view of the crucial role played by machining 
among all the material processing operations, a scientific 
understanding of various metal-cutting processes is vital for 
modern engineering technology and practice. 

A detailed analysis of the mechanics of metal cutting which 
sheds light on the underlying plastic defcormations, thermal 
phenomena and diffusion phenomena, is helpful in the prediction of 
cutting tool wear. 
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Over a past few decades, a large volume of experimental data 
has been collected on cutting tool wear, but a comprehensive 
theoretical understanding of this highly complex phenomenon is yet 
to emerge. In order to estimate the tool wear rate, one has to 
rely on empirical relations based on experimental measurements and 
material properties. 

The present work, in its right ernest, ’s a small effort to 
fill some existing gaps in the theory of cutting tool wear. A 
coupled analyses of heat transfer and diffusion throughout the 
whole region of the work piece, chip and tool has been attempted 
here, to provide a reasonably accurate estimate of the tool wear. 


1.2 REVIEW OF PREVIOUS WORKS 
1.2.1 MECHANICS OF CHIP FORMATION 

The earliest systematic study of metal cutting process, dates 
back to the end of 19th century, when pioneering researchers like 
Mallock, Tresca etc. attempted to throw some light on the 
mechanics of machining- However, a major contribution to this 
area which had an important influence on the practical development 
of machining the work by Taylor C19073. He was particularly 
concerned with toed wear and tool life and seems to be the first 
to recognize the influence of temperature on tool wear. He 
developed a relationship between tool life and the speed of 
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machining which is still in use today. 

After this work, research in the area of metal cutting 
commenced from the post war period. The first attempts towards 
developing an understanding of the mechanics of metal cutting were 
performed by analyzing the chip formation process. A pioneering 
work was carried out by Merchant [1944, 19453 on the so called 
classical "single shear plane model". This model assumes that the 
chip is formed due to plastic deformation occuinring along a single 
plane, known as the 'shear plane'. The most advantageous aspect 
of this model is that it enables the determination of the average 
yield shear stress and slip velocity along the shear plane, purely 
through geometric constructions. It also established beyond any 
doubt that metal cutting is basically a shear deformation process. 

Piispanen £19483 performed an analysis similar to that of 

Merchant £19453. He hzis considered the additional aspects of 
internal friction and strain hardening. He has shown that a high 
cutting speed decreases the friction angle on the face of the 
tool, by which the forces loading the machine tool steeply may be 
decreased. 

Drucker £19493 investigated two dimensional orthogonal 

cutting in detail, introducing several concepts and hypotheses 
which apply equally well to the most general cutting operation. 
Of these the most significant is the high rate of plastic 

deformation process. He has explained Ekstein's paradox and 
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correlated the influences of speed of cutting, depth of cut and 
rake angle by the consideration of dynamic plasticity. In 
addition, the effect of non-homogeneity of the material is 
analyzed and shown to account for the gradual change from the 
discontinuous to the continuous type chip as the cutting speed is 
increased. 

The shear plane model of Merchant [1944, 19453 and Piispanen 
[1948] was improved by Palmer & Oxley [19543 and Okushima & Hitomi 
[19573. These studies suggested that through the formation of 
chip occurs due to the flow of metal under shear and deformation 
is not limited to a single plane; it occurs in a narrow region 
approximated by two parallel planes. This model later came to be 
known as the thick shear zone model. In view of the small 
thickness of the shear zone in comparison to its length, this 
model assumes that the state of stress within the deformation 
region is uniform simple shear. 

Around the same period of the above mentioned theoretical 
models, efforts were going on to validate them with 
experimentation . Kececioglu [19583 was the first person to 
observe the process through photo-micro graph. The plastic zone 
where the chip formation takes place was photographed at short 
intervals using a quick stop device. By noting the deformation of 
the grain lx)undaries, the size and the shape of the plastic 
deformation zone were estimated. This study established that the 
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plastic region could be approximately represented by a thin 
parallel-sided zone. 

Later, a more effective way of using photo- micro graph 
technique was proposed by Stevenson & Oxley C19693. Assuming a 
thin parallel-sided zone, they observed printed grids on the 
material before and after deformation. From the streamlines of 
metal flow, the strains and strain rates were calculated. 

1.2.2 PREDICTION OF TEMPERATURE IN ORTHOGONAL CUTTING 

Hann [19513 calculated the shear plane temperature by 
assuming the shear plane to be a uniform band of heat source 
moving obliquely through an infinite work piece. 

Loewen & Shaw [19543 also assumed the shear plane to be a 
uniform band of heat source but considered that it moved over a 
semi-infinite work piece with the proportions of the shear plane 
heat entering the work piece and the chip determined by Blok's 
partition principle. 

Trigger & Chao [19513 used Blok's [19383 partition principle 
to calculate the average tool-chip interface temperatxire assuming 
that the heat dev^oped at the tool-chip interface was uniformly 
distributed. 

Weiner [19553 obtained a solution for the shear plane 
temperature distribution by assuming that the chip velocity was 
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perpendicular to the shear plane and that the heat conduction in 
the direction of motion of the work piece and chip could be 
neglected. 

Rapier [19543 calculated temperature distributions in the 
work piece, chip and tool which he treated as three separate 
systems. He also assumed a constant shear plane temperature and a 
plane uniform heat source at the tool-chip interface with all the 
heat flowing into the chip and none into the tool. 

Dutt & Brewer [1964] improved the analysis by treating the 
work piece, chip and tool as one system, but after making some 
approximations found that they were able to dispense with the tool 
region altogether. In this way they were able determine the 
proportions of shear plane heat entering the chip and tool. 

Chao Sc Trigger [19553 improved their earlier analytical 
solution for the interface temperature distribution by allowing 
the fraction of interface heat flowing into the tool to vary along 
the interface although still assuming a uniform heat source. 

The disadvantages of the above methods of temperature 
calculation result from the simplifications made in the shear 
plane model. This assumes a velocity discontinuity along the 
shear plane while in real situations, the transition from the work 
piece to chip velocity occurs gradually over a finite plastic 
2one. Also the velocity of the chip material adjacent to the 
tool-chip interface (i.e. viithin the secondary zone ) is less than 
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the chip velocity, resulting in the characteristic deformation 
observed in this region. 

It has been assumed in all the cases that the generated heat 
is unifomly confined within the shear plane and the tool-chip 
interface rather than being spread over the finite primary and 
secondary plastic zones. A further disadvantage is that, with the 
exception of the work of Dutt & Brewer [1964], the shear plane & 
interface temperatures were calculated independently assuming that ’ 
there is no interaction between the two heat sources when in fact 
some interaction would be expected. 

1.2.3 VISCOPLASTIC MODELS: DEFORMATIONS. TEMPERATURE PREDICTION 
USING FINITE ELEMENT METHOD 


The credit for the effective usage of the Finite Element 
Method (F.E.M) to viscoplastic analysis goes to Zienkiwicz et al~ 
[1981]. In 1978, these researchers applied FEM for the modelling 
of metal forming and extrusion problems. In metal cutting, the 
FEM technique was employed by Tay & Stevenson [1974, 19763 for 
predicting the temperature field for orthogonal machining. They 
used hyperbolic streamlines in the primary deformation zone for 
calculating the velocity field & strain rates. They also derived 
semi-empirical relations for computing heat generation. However, 
these studies do not provide a coupled analysis of plastic 
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deformation and thermal processes during metal cutting. 

Stevenson, Wright & Chow [19833 modified the work of Tay, 
Stevenson & Davis Cl 9763 in calculating the temperature 
distributions in the chip and tool during metal cutting for an 
extended range of shear angles & contact lengths. They observed 
good agreement between their FEM results and earlier experimental 
measurements. 

Murarka et al. [19793 predicted the temperature field 
variation during metal cutting and verified the average 
temperature rise by experimental measurements. These authors 
obtained the temperature distribution near the tool-tip for a wide 
range of cutting conditions. 

Balaji et al. [19863 applied the FEM technique to evaluate 
the temperature field within the work material and the tool, for 
coated carbide tools. The effects of protective coating upon the 
tool-chip interface temperature, cind heat transfer characteristics 
have been investigated in their work. They have also provided the 
experimental verifications of the results obtained. 

Strenkowski & Carroll [19853 used FEM to predict the chip 
geometry, plastic deformation and the residual stresses in the 
work piece. They have presented an adiabatic heating model to 
stimulate the heat generation effects due to plastic work. But 
they did not consider the coupled temperature field in the tool 

Strenkowski & Moon Cl 9903 modified the work of Strenkowski St 
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Carroll [1985] by considering the temperature distribution in the 
work piece, chip and tool without the need of empirical cutting 
data. 

Sarma El 990] calculated the nature of temperature 
distributions in the work piece and the chip during machining. He 
assumed that heat is generated due to plastic deformation at the 
primary zone and due to plastic deformation as well as frictional 
rubbings at the secondary zone. He considered constant values for 
material properties such as thermal conductivity, specific heat 
and yield stress. 

Despande [1992] calculated the temperature distributions in 
the work piece and chip during machining using the similar 
concepts as Sarma [1990]. But he considered the variation of 
thermal conductivity, specific heat and yield stress with 
temperature in calculating the temperature field- 

In the most recent times, Komvopoulos et al. [1991] have 
analyzed the deformation processes in metal cutting using 
elasto-plastic approach. They did not consid^ the thermal 
aspects of the cutting process. 


1.2.4 CUTTING TOOL mAR 


Numerous experiments have been conducted on tool wear during 
the past few decades and much is still to be done. No general 
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laws of tool wear have been deduced so far. Since metallic wear 
is, in general, a complex phenomenon involving chemical, physical 
and mechanical processes, it is clear that a theoretical treatment 
of the problem of tool wear at present must be accompanied by 
simplifications. 

Dawihl [1941], after a series of investigations on crater 
wear of cemented-carbide tools, concluded that mutual diffusion of 
materials at the tool-chip interface and the propensity to welding 
of certain constituents in sintered carbide with steel were of 
direct concern. He pointed out the importance of temperature at 
the tool-chip contact. 

Trent [1952] proposed a new theory to account for the crater 
wear of cemented-carbide tools. He postulated that fused layers 
of alloy were formed between the steel and the free 
tungsten-carbide in the cutting tool material at the tool-chip 
interface. Wear was the consequence of the transportations of the 
fused alloy by the moving chip. The extent of cratering was 
controlled by the distribution of temperature near the cutting 
edge- 

Trigger & Chao [1956] also found that wear at the top surface 
of the tool is essentially of the transfer type and formation of 
crater wear is strongly temperature dependent. 

Cook & Nayak [1966] presented an empirical relation to show 
that the geometry of the tool wear is strongly controlled by the 
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chip curl. They presented simple theories to calculate the rate 
of tool wear based on diffusion of tool material into the work 
piece. 

Rubenstein [1975] deduced an analysis of continuous flank 
face wear which is attributed only to adhesion between tool and 
work piece material. His analysis is based on the assumptions 
that flank face wear occurs exclusively by the adhesion mechanism 
and 'weld' formation occurs by a diffusion process which is 
temperatvire and time dependent. 

Kannatey-Asibu [1985] derived a two dimensional 
transport-diffusion equation for in metal cutting, by an extension 
of the Fick's second law of diffusion. The analysis considers the 
dynamics of the cutting process and thus offers a comprehensive 
representation of the diffusing species in the work piece and the 
chip. 

1.3 OBJECTIVES AND SCOPE OF THE PRESENT WORK 

In the present study, an attempt has been made to fill the 
gap in the analysis of cutting tool wear by solving the coupled 
heat transfer and diffusion equations- 

It is well established from previous experiments that the 
deformation is restricted to a very small region around the 
tool-tip. Therefore, a solution domain extending upto a few 
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millimeters from the tool-tip into the workpiece and the chip has 
been considered. Beyond this domain, work material is assumed to 
be perfectly rigid. The tool is also taken to be perfectly rigid 
and sharp, with no built-up edge. The tool cutting edge is 
assumed to be orthogonal to the cutting direction. Heat 

generation during machining is restricted to the primary and the 
secondary zones only. The heat generation in the primary zone is 
due to plastic deformation only whereas in the secondary zone, it 
is due to both plastic deformation and interfacial friction. 

The changes in material properties such as the specific heat, 
thermal conductivity and the diffusion coefficient with 
temperature have been taken into account, since metal cutting 
involves large temperature variations. Based on the temperature 
field within the whole domain, carbon concentration has been 
calculated throughout the region considering two-dimensional 
species-diffusion eguation. The yield stress for the tool 
material has been assumed as a function of both the temperature 
and the carbon concentration. The actual shear stress which is 
exerted by the chip on the tool surface, varies with distance from 
the tool-tip. From the tool-tip to the 'sticking ' portion, it is 
constant and beyond that it decreses linearly to zero at the point 
where the chip leaves the surface. In the present work, it has 
been assumed that if actual s\u:face stress is nK>re than the 
calculated yield stress - erosion occurs resulting in 'crater' 
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wear on the rake face. Prediction of crater wear profile has been 
done on the basis of this assumption for two cases of cutting 
conditions. 

A special software for FEM grid generation has been 
developed, since grid generation for the complex geometry 
encountered in the present study is quite a difficult task. Some 
of the advanced concepts such as the hyperbolic streamline model { 
lay et al. [1974] ) and the trans-finite interpolation have been 
used in the grid generation algorithm. 

1.4 ORGANISATION OF THE THESIS 


Chapter 1 of the thesis provides the introduction of the 
problem and a review of literature relevant to the basic 
mechanisms involved in orthogonal metal cutting which influence 
the temperature field around the tip and the cutting tool wear. 
The theoretical formulation of the problem is presented in Chapter 
2. Chapter 3 details the finite element analysis of the heat 
transfer and diffusion equations in the vicinity of the cutting 
edge. Chapter 4 presents the numerical results. Comparison of 
the resxalts obtained with those of the earlier researchers has 
been included in this chapter. Chapter 5 presents the conclusions 
cind the scope of further research. 



CHAPTER 2 


THEORETICAL FORMULATION 


Metal cutting involves very complex deformation phenomena. 
Extremely localized asymmetric deformation occurs at exceedingly 
high strain rates, in a minute domain of few millimeters size- 
The metal cutting action itself is a consequence of high 
compression followed by shearing of the metal. Consequently, any 
theoretical treatment of the metal cutting process is bound to 
involve certain relaxations and simplifications. 

A two dimensional orthogonal cutting condition has been 
considered for the analysis. In fact, orthogonal cutting is a 
particular case of oblique cutting, with cutting edge of the tool 
perpendicular to the cutting velocity vector. The tool is assumed 
to be sharp, with no built up edge. The work material is taken to 
be ductile, resulting in the formation of continuous chips. The 
tool and the portion of the work material outside the solution 
domain on the other hand, are taken to be perfectly rigid without 
any deformation. 


2.1 PLASTIC DEFORMATIOli & CHIP FORMATION IN ORTHOGONAL CUTTING 
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In orthogonal cutting, the fact that the width of the chip is 
large compared to its thickness, renders the problem two 
dimensional. Although the deformation in metal cutting is similar 
to those of the metal forming processes, a distinguishing feature 
is that the deformed material, is separated from the parent work 
piece in the form of a chip. Hence, it is of primary importance 
to analyze the chip formation process, in any theoretical analysis 
of metal cutting. 

The uncut material in motion is first interrupted by the 
cutting tool, causing severe compressive stresses in the material. 
After sufficient compression, when the resulting shear stresses in 
the material attain the yield value, plastic flow of the material 
occurs in the direction of shear force leading to the formation of 
the chip. The outward shearing movement of each successive 
deforming element is arrested by work hardening and the movement 
is transferred to the next element. 

It is wen accepted now that the predominant mechanism of 
metal cutting is that by shearing process occurring in the zone 
called ’shear zone'. A detailed study of metal cutting reveals 
that two such 2 ones exist in the domain. 

(i) Primary shear deformation zone (PSDZ) 

(ii) Secondary shear deformation zone (SSDZ) 

These two zones have been shown in Fig. (2.1). The primary 
shear deformation zone has been observed to be wedge shaped and 
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Fig. 2.1 ; Shear Deformation Zones 
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its size depends upon cutting velocity, it has been reported that 
an increase in cutting velocity causes the thickness of PSDZ to 
reduce. The thickness or mean width of PSDZ is an important 
quantity which gives an idea of the rate of deformation and the 
amount of heat generated. The narrower the region, the greater 
the rate of deformation. 

The secondary shear deformation lies adjacent to the contact 
patch between the chip and the tool. Additional deformation due 
to friction takes place here. Each layer in SSDZ moves at a 
different speed because of shearing in a direction parallel to the 
rake face of the tool. An interesting point to note is that the 
highest temperature occurs on the interface between the chip and 
the tool. The temperature in SSDZ is higher than that in PSDZ 
mainly due to the fact that the metal flow velocity is very low in 
SSDZ and also that the material which enters SSDZ has already been 
heated through PSDZ. Some of the researchers (Jain et al. [1987]) 
have attempted to estimate the thickness of these zones 
theoretically and experiraentally- 

2.2 GOVERNING EQUATIONS 

2.2.1 TEMPERATURE FIELD IN THE CUTTING DOMAIN 

In the present analysis, temperature fields in the plastic 
deformation zones of the work piece, chip and the tool during 
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machining are analyzed. It is obvious that in order to get a real 
picture of the temperature distribution in the whole region around 
the cutting edge, heat generation must be calculated carefully. 
It has been assumed that all of the work done is converted into 
heat, which seems to be quite reasonable because the material 
density does not change considerably and also, work hardening is 
not severe due to the high temperature conditions. 

The steady state heat balance equation in the regions 
surrounding the cutting edge can be written in a general form 








ax 


+ V 



+ Q,= 0, 


( 2 . 1 ) 


where T is the temperature field, u and v are the components of 
the velocity field, k, p, c are the thermal conductivity, density 

p 

and specific heat respectively, and is the rate of heat 
generation per unit volume. The subcript i is used to denote 
different regions around the tool tip such as the fully formed 
chip, uncut material, machined region, PSDZ, SSDZ and the tool. 
It is evident that in order to solve equation (2.1), it is 
necessary to find out the materM properties (k, p, c^), the 
velocity field (u, v) and the rate of heat generation Q in each 
zone. The procedures followed to calculate each of the above 
quantities are detailed in the following sections. 
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2.2.2 CALCULATION OF VELOCITY FIELDS 

For calculating the velocity jEield at different locations 
some approximations have been made. Within the uncut metal and 
fully formed chip, the velocity field is taken to be uniform given 
by V and V respectively, where V = { irD N ) / 1000 and V 

Chip 1 chip 

Vr. In the above expression, r is the chip thickness ratio given 
by r= t/t 

C 

In the primary zone (I) and the secondary zone (IT) (Fig. 
(2.1)) where plastic deformation takes place, velocity varies from 
location to location. For the primary zone the velocity at any 
location is obtained by linearly interpolating between the values 
V and The primary zone is divided into two layers of equal 

thickness for the sake of velocity calculation. The elements in 
the first row (I^) have velocity components of the form : 

u = -(V . . cosa ) / 3. 

chip 

I 

and,v = -(2.V +V sina) 3. j (2.2 a, b) 

chip 

The elements in the second row (I^) have velocity components 
given by : 



and,v 


sina) / 3 


(2.3 a, b) 
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The velocity components in the fully formed chip are t 


and. 


c 

il 

1 

< 

cosa 

chip 

> 

V = -V _ 

sina 

chip 



(2.4 a, b) 


The components of the metal flow velocity in the secondary 
zone depend upon the sliding v^ocity of the chip at the 
interface. Tay et al. [19763 found that the sliding velocity of 
the chip at the interface is given approximately by the equation - 

V. = 4 (V ) / (1 + (8 X) /h), (2.5) 

mi o chip 

where 'x' is the distance from the cutting edge measured along the 
interface and 'h' is the total contact length. 

For calculating the velocity field in the secondary zone, 
linear interpolation is used between the chip velocity and 

the interfacial velocity (V. ) similar to that in the primary 

mt 

zone. In the tool region, the velocity components u and v are 
taken to be zero. 


2.2.3 HEAT GENERATION DURING MACHINING 

The rate of heat generation per unit volume is assumed to be 
zero everywh^e except in the primary and secondary deformation 


zones. 
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Heat generation in the primary shear deformation zone is due 
to plastic work only. The rate of heat generation (per unit 
volume) at any point in the primary zone is given by 

Q = (fs vs) / volurael, (2.6) 

pz 

where fs, vs and volumel are the shear force, shear velocity and 
volume of the primary deformation zone respectively. Here, shear 
force can be calculated using Merchant's circle diagram (see, Fig. 
(2.2)), using the expressions : 

fs = fh cos (p - fv sin 
and, vs = V cos a / cosi<P - a), 

where V, 4> and a are the cutting velocity, shear angle and rake 
angle respectively. Volume of the primary deformation zone is 
given as the product of the 'width of cut’ and the area of PSDZ. 

In the secondary shear deformation .zone, there are two types 
of heat generation. One is due to the secondary plaistic work and 
the other is due to the interfacdal friction. Upto to the point 
of 'sticking friction' both types of heat generation are present, 
whereas beyond the point of sticking friction only frictional heat 
generation is present. Shear stress is constant in the sticking 
friction point and after that it varies linearly with rest of the 
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contact length and becoming to 2 ero at the point where the chip 
leaves the tool rake surface (see, Fig. (2.3)). 

Heat generation due to secondary plastic work is given by 

Q = f2a (V — V ) / deita2 (1- x /h), (2.7a) 

cn 1 p i n t 

where f2a, h and delta2 are the shear stress, the contact length 
and the thickness of the secondary deformation zone respectively. 
Also, X is the distance from the cutting edge along the contact 
length. Here, the SSDZ is assumed to be of triangular in shape 
{Tay et ai.[1976]), consisting of one layer of elements adjacent 
to the rake face within the chip. Its maximum thickness delta2 
has been taken as ; 

delta2 = (0.154 t) / sin (2.7b) 

Heat generation due to frictional work been calculated 

as follows 


= f2b V. . / delta2 (1- x /h), 

2b int 


( 2 . 8 ) 


where f2b and V are the shear stress and the interfacial 

int 

velocity (see, equation (2.5)) respectively. 
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Fig. 2.3 : Variation of Shear Stress on the Rake Face of the Tool 
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2.2.4 EVALUATION OF THERMAL PROPERTIES 

A careful study of the governing equations reveal that the 
material properties such as thermal conductivity and specific heat 
affect the temperature solution to a considerable extent. All the 
above mentioned properties vary with temperature. In particular, 
at high temperatures, the variations become quite significant. 

In the present analysis, the solution domain is a minute 
region surrounding the tool tip. As already mentioned elsewhere, 
the temperatures in this region are quite high due to the rate of 
plastic deformation. Thus, to have a realistic picture of the 
results, the variations in these properties at higher temperatures 
are incorporated in the analysis, where material data is 
available. 

The variation of the thermal conductivity and the specific 
heat are curve fitted and property correlations have been derived. 
For thermal conductivity and specific heat , piece wise equations 
of a straight line of the type y = rox + c have been fitted for 
different temperature ranges. The data used in the present study 
is shown in Figs.' ((2.4), (2.5) and (2.6)). 
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Fig 2.4 : Variation of Thermal Conductivity with Tempearture in 

Workpiece (Mild Steel) 
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Fig2.5 : Variation OF Thermal ConductivityvithTempearture IN 

Tool (High Speed Steel) 
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2.3 ANALYSIS OF TOOL WEAR 

A cutting tool may fail because of any one of the mechanisms 
such as the process of softening due to plastic deformation, 
brittle fracture due to mechanical loads, shocks and gradual wear. 

The failure due to plastic deformation of the cutting edge 
takes place when the hardness of the cutting edge is not 
sufficient to maintain the original shape of the tool. The 
stresses developed because of the resistance of the work material 
at the prevailing temperature and strain rate cause rounding of 
the cutting edge. According to Bhattacharya & Sen E19841, it has 
been shown that if the ratio between the hardness of the tool and 
the chip material exceeds 1-5, such failures are not expected. 

Failures by brittle fracture are due to excessive pressure 
and load or shocks such are present in intermittent cutting, quick 
freezing of the cutting action, cutting under chatter condition, 
etc. However, brittle fractures are avoidable by proper selection 
of cutting conditions (feed, depth of cut, etc.) or by 
strengthening the cutting edge with a land. 

But still, cutting tools continue to fail by a process of 
wear which is due to the interactions between the chip and the 
tool or between the work piece and the tool. After the cutting 
has progressed for some time, wear which appears on the tool face 
because of tool-chip interaction, forms a cavity known as 
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'crater'. It begins at a certain distance from the cutting edge. 
On the other hand, wear which takes place due to the rubbing of 
the tool with the work piece on the flank face of the tool is 
known as 'flank wear*. 

2.3.1 MECHANISM OF WEARING PROCESS 


Wear is not a process involving one unique mechanism — there 
are several mechanisms which are known to be operative, depending 
on the physical conditions existing. Four main wear mechanisms 
have been identified as follows : 

(1) Adhesion 

(2) Abrasion 
(31 Diffusion 
(4) Fatigue 

(1) Adhesion Wear : Due to the existence of a very high 
temperature region on the chip-tool interface, formation of 
temporarily welded junctions between the tool and the chip and the 
subsequent destruction of these occurs continuously. When the 
destruction is by shearing below the interface, a wear particle is 
transferred. This mechanism is termed as the adhesion wecur. 

(2) Abrasion Wear : The abrasion process involves cutting at the 
microscopic level and as such, it depends on the hardness, the 
elastic properties and the geometry of the two mating surfaces. 
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It has been pointed out that the larger the amount of elastic 
deformation a surface can sustain, the greater will be its 
resistance to abrasion. 

(3) Diffusion Wear : Holm proposed diffusion wear as a process of 
atomic transfer at contacting asperities which eventually causes 
material softening leading to wear. It is a highly 
temperature-dependent phenomenon, i.e., a direct function of the 
rubbing speed. However, the amount of material transferred by 
diffusion is also dependent on the time of contact of the mating 
surfaces, which is an inverse function of speed. 

(4) Fatigue Vfear : This is described as the process of change in 
stress in the material below a surface as an asperity from a 
rubbing surface passes over it. Each asperity is associated with 
a wave of deformation. At some distance ahead of the asperity the 
underlying material is compressed, but behind the asperity, 
tensile stress elongates the material. This chamge in sign of the 
stress as an asperity passes a given point causes fatigue failure 

of the material below the surface. 

In the present study, tool wear caused by the diffusion 
mechanism alone has been studied. The other mechanisms have been 
neglected for the sake of simplicity. 


2-3.2 DIFFUSION WEAR PHENOMENA 
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The fact that the workpiece and tool materials are different 
from each other causes diffusion of carbon and other alloying 
elements from one material to the other. This mass transfer 
occurs in both the directions. Close investigation near the 
tool-chip interface has revealed that : 

(1) There is constant migration of carbon from tool to chip. 
Since the carbide phase is responsible for tool hardness, carbon 
diffusion reduces the tool strength. 

(2) There are thermal gradients in both the tool and the 
chip. However, because the thermal diffusivity is many orders 
larger than the mass diffusivity, transient variation . of mass 
concentration needs to be considered while the temperature. field 
can be treated as steady. 

(3) The strength or hardness of materials is dependent on 
gross carbon diffusion and on the average, its value is higher in 
tool than in workpiece. 

Various factors which influence the diflEusion wear are : 

(A) GEOSS TOOL HARDNESS ; If at the operating temperature, the 
tool is soft enough that the interfacial stresses cause gross 
deformation, the tool geometry will clearly deteriorate. This 
imposes a hardness criterion on the tool material which may be 
separate cind distinct from a wear criterion. In steel tools, this 
mechanism can Inhibit cutting due to diffusion-controlled carbon 
depletion of a relatively deep surface layer or to 'annealing in 
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which the carbon diffusion permits a change in crystal structure. 

(B) GROSS DIFFUSION OF TOOL INTO THE WORK PIECE : In this 
mechanism, the tool matrix material, or a major strengthening 
component, is dissolved by the chip & work surface as they flow 
past the tool. The mechanism is somewhat analogous to the 
dissolution of a salt crystal in a moving stream of water. This 
effect could be present in carbide tools with carbon diffusing 
from the WC into the work. In cast alloys and ceramics, this 
could be predominant. In steel tool, it is possible that the iron 
could diffuse into the work, if the work material is non-ferrous. 

(C) DIFFUSION OF MINOR TOOL CONSTITUENTS WTO THE WORK : Here, a 
minor (but strengthwise important) component of the tool material 
diffuses into the work, leaving a thin surface layer with degraded 
physical properties. In tungsten carbide tools, cobalt which 
serves as the binding material could diffuse away, leaving the 
relatively fragile WC structure to break away. In ste^ tools, 
the carbon could diffuse away, forcing an alteration of crystal 
structure in a very thin surface layer- 

(D) DIFFUSION OF WORK MATERIAL COMPONENTS INTO THE TOOL t Here, 
either the work material or an additive element may diffuse into a 
superficial tool layer, altering its physical properties. It is 
well known, for instance, that many liquid metals (gallium, 
mercury, lead & so on) produce very strong embrittling effects on 
other materials. It is conceivable that a work additive such as 
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lead could diffuse into the tool surface, rendering a surface 
layer quite brittle. 


2.3.3 TRANSPORT DIFFUSION EQUATION IN METAL CUTTING 


The cutting process is two dimensional for orthogonal 
cutting, involving bulk motion of material in the direction of 
diffusion as well as in a direction perpendicular to it. Taking 
into consideration both the process of diffusion and the bulk 
motion, the equation governing the transient variation of the 
diffusing species can be written in the form 


dc §2 

at ~ axr 



(2.9) 


where, c is the concentration and u, v are the components of 
cutting velocity- These velocities are calculated in a similar 
way as was done for the evaluation of the temperature-field during 
machining. 

For solid media, the diffusion coefficient (D) varies with 
temperature in a form given by 

0 = D„exp 


where, is the pre-exponential constant, q is the activation 
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energy for diffusion and R is the absolute gas-constant. During 
computation, the values of = 0.7 * 10”^ /sec and q = 157 
KJ/mol were used for both HSS and carbide tools. 

It is possible to obtain the concentration field of any 
diffusing species, by solving equations (2.9) and (2.10). It 
would then be of interest to predict the deterioration in tool 
strength based on the temperature and concentration fields, and 
analyse its effects on tool wear. 

2.4 SELECTION OF THE SOLUTION DOMAIN FOR ANALYSIS 

It is well established from previous experiments that the 
deformation is restricted to a very small region around the tool 
tip. Therefore, a solution domain extending up to a few 
millimeters from the- tool tip into the work piece, tool and the 
chip has been considered. Beyond this domain, the material is 
assumed to be perfectly rigid. The tool is also taken to be 
perfectly rigid and sharp with no built up edge. The tool edge is 
assumed to be orthogonal to the cutting direction. Further, for 
generating the domain & the finite element mesh easily, the 
boundary of the domain has been divided into ten surfaces, as 
shown in Fig. (2.7). 

Surface I is the top portion of the tool extending from the 
tool shank to the end of contact length. 
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Surface n is the portion of the tool where it is clamped. 

Surface m is the bottom portion of the the cutting tool. 

Surface IV is the flank face of the cutting tool. Since the 
tool edge is sufficiently sharp and the clearance angle is 
adequate, there is no contact of the flank face with the machined 
portion of the work. 

Surface V constitutes the machined surface and is exposed to 
the atmosphere due to the clearance angle of the tool. 

Material leaves the plastic zone at surface VI. 

Surface vn is selected in a plane parallel to the cutting 
velocity and restricts the boundary of the plastic zone below the 
depth of cut. 

Surface vm is the entry section where material enters the 
plastic solution domain and is considered to be plane, 
perpendicular to the cutting velocity- 

The free surface of the chip, where the complex chip curling 
effect with serrations is observed, has been uKDdeled using the 
hyperbolic stream line concept as used by Tay et al. [19763. This 
forms the surface IX. 

The chip is considered to be bounded by a plane parallel 
surface beyond the contact point, which is represented by surface 
X. 


2.5 BOUNDARY CONDITIONS 
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2.5.1 TEMPERATLfRE FIELD 


Referring to Fig. (2.7), on surfaces (2), (3), (6), (7), (8) 
the temperature is prescribed as that of the ambient. The other 
surfaces ( (1), (4), (5), (9), (10) ) have convective boundary 
condition which is expressed as 


-k 


an 


h (T - T ) 

amb 


r 


where h'" is the film heat transfer coefficient and n is the normal 
coordinate for the surface under consideration. The value of h^ 
has been chosen as 50 w/mK while T ^ is set equal to 293 K. 

amb 


2.5.2 CONCENTRATION FIELD 


Surfaces (2), (3) have a concentration of carbon is equal to 
0.85 % and surfaces (7) and (8) have carbon concentration is equal 
to zero (see. Fig. (2.7)). For other surfaces, the mass transfer 
boundary condition has been taken as 

which implies that these surfaces are impervious (zero flux). 


2.6 INITIAL CONDITIONS 
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The initial condition of temperature field for the tool, 
workpiece and the chip has been taken as that of ambient 
temperature. 

Initial condition of concentration field for the tool has 
been set equal to 0.85 % and for the workpiece and the chip it is 
taken as 0 %. 

2.7 WEAR PREDICTION 

The governing equations for the temperature and concentration 
fields (equation (2.1) and (2.9)) have been solved numerically 
using the finite element method, in the present work. From the 
predicted temperature and concentation fields, changes in strength 
(yield stress) have been estimated for the tool material at each 
location. Eventually, based on these data, the expected tool wear 
profiles have been predicted. The procedure for estimating the 
wear profile is discussed below. 

In the present work, it has been assumed that Ccirbon is the 
only diffusing component from the tool to the workpiece. Although 
this is too simplified a picture for realistic predictions, all 
the important diffusing components can be taken into consideration 
in a similar fashion. For HSS tools, the influence of carbon 
concentration on tool hardness (standard Rockwell C) has been 
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provided in graphical form by Clark et aZ.C1987] { see Table (2.1) 

). It is well known that the hardness depends on the temperature 
also. The data for the effect of temperature on the hardness has 
been obtained from Hoyle [19883 ( see Table (2.2) ), for type 
HSS tools. Finally, the hardness at any temperature and carbon 
concentration has been obtained in the form 

H(t,c) = h{T) * f{c), (2.11) 

where f(c) is a function whose value is between 0 and 1. For C > 
0.60%, f{c) = 1.0. The function h(T) gives the variation in the 
Rockwell C hardness with respect to the tepmerature at a standard 
carbon concentration of 0.85%. Using stcindard tables for steels, 

{ Clark et aZ. [19873, see Table (2.3) ) the Rockwell C hardness 
has been- converted into equivalent tensile strength 
FinaUy, the yield strenght in shear (u^^) has been estimated as ; 

a = a //3 (2.12) 

The variation of the actual shear stress on the tool face was 
presented in Fig. (2.3). By comparing the yield stress with the 
actual surface stress, it is assumed that wear is likely to occur 
whenever yield stress is less than he actual stress, at any 
location, on the chip-tod interface. In this manner, the length 
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of the crater on the tool surface is approximately predicted. 

The numerical solution procedure and the results obtained for 
some typical machining conditions are presented in the subsequent 
chapters. 


Carbon concentration 

{ % c) 

Factor, f{c) 

0.1 

0.6 

0.2 

0.6293 

0.3 

0.7846 

0.4 

0.8769 

0.45 

0.923 

0.5 

0.9615 

0.55 

0.9846 

0.6 

1-0 

0.7 

1.0 

0.75 

1.0 

0.8 

1.0 

0.85 

1.0 


■ Table 2.1 : Effect of Cahbon Concentkation on Hardness 
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Temperature (K) 

Hardness, h{T) 

321.0 

65.0 

368.0 

63.8 

588.0 

62.0 

698.0 

60.0 

813.0 

57.0 

923.0 

. 47.0 


SMPERATURE ON HARDNESS 
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Table 2.3 : Conversion of Hardness to Tensile Strength 



CHAPTERS 


FINITE ELENENT ANALYSIS 


3.1 INTSODUCTION 


The finite element technique, whose birth and boom in the 
1960s was due to the application of digital computers in 
structural analysis, has spread to a variety of engineering & 
physical science disciplines in the last two decades. In many 
engineering problems, it is not always possible to obtain a closed 
form exact solution. Thus, for such difficult problems, it is 
possible to get the best approximate solution (i.e. very close to 
the exact one), using a powerful numerical technique such as the 
finite element method. 

The basic concept of finite element method is that of 
converting a continuous field problem into an approximately 
equivalent discretized problem- A finite number of points are 
identified in the problem domain and the values of the field 
variable to be solved { and its appropriate derivatives ) at these 
■ points are selected as unknown variables. The points are called 
as nodes. The domain of the function is represented approximately 
by a finite collection of sub - domains called elements, which are 
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connected together. The solution variable is approximated locally 
within each element by continuous functions that are uniquely 
described in terms of nodal point values associated with the 
particular element. The path to the solution of a finite element 
problem consists of essentially five specific steps : 

(i) integral formulation of the problem, 

(ii) discretization & definition of the elements, 

(iii) establishment of the nodal equations for one element, 

(iv) assemblage of elemental contributions to nodal equations 
and incorporating boundary conditions, and 

(v) the numerical solution of the global matrix equations. 

The formation of nodal equations can be done using any one of 

the procedures such as the direct approach, variational method, 
method of weighted residuals, and the energy balance approach. 

In the present analysis, the Galer kin’s weighted residual 
method has been used because of its versatility. 

3.2 APPLICATION OF FEM 
3.2.1 GEID GENERATION 

The first step in the finite element method is the selection 
of solution domain and its discretization. The choice of the 
solution domain boundaries for the metal cutting problem was 
discussed in the previous chapter. Grid generation for such a 
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domain is quite a tedious task. This is in view of the that the 
variations in depth of cut, chip thickness, rake angle or contact 
length, lead to totally different shapes for the solution domain. 
Thus, it is obvious that unless one has an efficient way of 
generating grid, the analysis can not be done conveniently. Grid 
generation, therefore, forms one of the important aspects of the 
present analysis. 

Coordinate data and nodal connectivity of elements for each 
2one are calculated separately first and global element and nodal 
numbers for all 2 ones are assigned later. To begin with, the 
coordinates for all the boundary nodes are prescribed suitably and 
then using a special type of interpolation technique (known as the 
trans-finite interpolation) coordinate data of interior points 
have been generated. This technique is based on transforming an 
arbitrary 4-sided domain into a square-domain, in terms of 
body-fitting coordinates (C, n) as shown in Fig.(3.1). 

The procedure for generating the X and Y coordinate values 
for a typical point H is illustrated in the figure. The domain 
ABCD is mapped onto a unit square in the transformed (C, h) plane. 
If 'H' is the int^section point of the grid lines PQ & RS (which 
are t) » constant, and ^ = constant lines respectively), by 
bi-directional interpolation, the X & Y co-ordinate of H can be 
obtained in terms of the X & Y co-ordinates of the boundary points 
P, Q, R, & S. For instance, if C (1-0 are the distances of 
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Fig 3.1 : Body-fitting Coordinate Transformation 
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H from P and Q respectively in the transformed domain, then for 
one-dimensional interpolation : 

Xy= (1-D Xp + C (3.1a) 

Y^= {1-C)_ Yp + C (3.1b) 

However, for two directional interpolation, the 
interpolations must be applied in both C and h directions suitably 

then obtained in terms of the coordinates of the 
boundary points P, Q, R & S. Similarly, for obtaining the 
coordinates of aU the interior points, the appropriate boundary 
points on the 4 sides are selected and by trans-finite 
interpolation, the interior coordinates are generated. A typical 
mesh thus generated, is shown in Fig. (3.2). 

In the present analysis, a grid with 678 nodes and 604 
elements has been used, after the observation that a smaEer 
number of elements is not sufficient to accurately represent the 
sharp variations in the temperature and concentration fields. 
Provision has been kept in the program, to change the number of 
elements and the size of the elements as desired, to overcome such 
problems. 

3.2.2 derivation OF FINITE ELEMENT EQUATIONS 
(i) ENERGY EQUATION 
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FOUR NODED ELEMENT SHOWING THE LOCAL COORDINATE SYSTEM 

Fig 3.2 : Sample Grid Produced by Mesh Generator 
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The two dimensional heat transfer equation was presented in 
the previous chapter as ; 



-pc. 


9T , ar 

u a— + V a- 

i ( ax ayj 


+ Q =0 


(3.2) 


In order to obtain a discreti 2 ed solution, the temperature at 
any point within the region can be approximated as : 


{N} 


where IN) is the shape function array which depends upon the 
coordinates and is the nodal temperature vector. 

The residue at a point can be l^^ritten as. 




axl 



^ a f. ar^®’! , f 

raj J 

ayi^ay J ^ 


aT^-'^ 


ax 


+ V 


dr 

ay 


U)^ 


+ Q (3.3) 


Now by Galerkin’s weighted residual technique {Reddy [19863), the 
equation for the i^** node is obtained as 


IS N.R^^^dx dy = 0, 


(3.4) 


where N^ is the shape function corresponding to node i. Now we 


Coin write *. 
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IS N. 


a f 1 + a r „aT'-> 1 _ 

9x1^ 9x J 9y|^ 3y j 


P c. 


u 


9x 


(e) 


+ V 


at 

9y 


Ce) 


] 


+ Q dx dy = 0 


(3.5) 


Applying Green's theorem, the above equation reduces to the form 


r ^ aN . 9 n dN . 9N . <1 

"K ftT 


^ ^ I + P C N 

dy dy J P 1 


aN. dN. 'ii , . 
ai' + ^ a?' JJ{T}'""^dxdy 


^ aN . oiM . \ 

IS N,Q dx dy + # H^k[ n^+ gj' j 


aN. 


^ ri © ) *1 

{T} dB = 


= # N^k + J/ N. Q dx dy (3.6) 

The second term on the RHS is the heat generation 
contribution and the first term is from heat transfer boundary 
conditions. 

For evaluating the integral expressions in equation (3.6), 
Gauss-Legendre numerical quadrature is used. The numerical 
integration can be performed easily by introducing local 
coordinates CCr H) as shovsin in Fig. (3-1). In terms of these 
coordinates, any quadrilateral element is mapped onto a square 
with C and q variations over the range, ±1. In the present work, 
4 - noded quadrilateral elements have been used for which the 
shape functions are given by Fta 
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N,= J [l + (l + >7,n), (3.7) 

and, i == 1, 2, 3, 4 within each element. 

(ii) SPECIES-DIFFUSION EQUATION 

The diffusion equation is given by *. 



In a similar manner as done in the case of energy equation, 
the residue equation is obtained as : 


_ Ce) dc 
^ ^ dt 


(e) 



{3.9) 


where, 

{N(x,y)} {C(t)}^"®^ (3.10) 


The weighted residue equations are given by : 



(3.11) 
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Applying Green's theorem, the final form of the residue 
equations are obtained as : 


If N^N^{C}‘”®’dx dy + // D 


dNj dN. dNj ^ 

^ dx dx ^ dy dy j 


dx dy + 


f dN . dN 

^ dy{C}'“'’'= ^ N 


(ne) 




3r* 

^ dB 


(3.12) 


The transient term { first term on the LHS ) in equation 
(3.12) is handled by finite difference approach. There are 
essentially three types of discretization schemes for marching in 
time. These are : 

(i) Implicit scheme 

(ii) Explicit scheme 

(ih) Crank-Nicolson scheme 

In the implicit procedure, the time derivative C is evaluated 
at t ,, while marching from the n^** time level (t ) to the 
(n+l)^^ time lev^ (t ,). For the explicit scheme, C is 

n+1 

evaluated at t , while for the Crank - Nicolson it is evaluated 

n 

at t^^j. . A generalized implicit / explicit time discretization 
can be employed by using the expression : 


C = (1 - e) {Cl. + e {C}, 




(3.13) 
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where e is a fraction between 0 to 1. 

The transient discretized species finally obtained as : 



[MASS] + e ESTIFFNESSlIcc}^"!^ 

j n+1 


= At {FI -1- [MASSICC}^"®^ 

n 


- (l-e) At [STIFFNESS3{C}^"®\ 

n 


(3.14) 


where. 


[MASS] = IS N.N^dx dy 



r 8fli dN. 

1 L 

( 6x dx 

dN, 

aH ^ 

] ■ax ay 

[STIFFNESS] = ff D 

dy 

f 

dN 

dN. 

55^ 


+ // N. U 

+ V 

dx 

dx dy 


and, {F} = # N.l^ dB 


During computations the implicit scheme has been used and 
hence e =1. The time step (At) has been chosen in the range (5 to 
20) seconds. 

3.3 BOUNDARY CONDITIONS 

Modelling of boundary conditions is as important as the 
modelling of differential equations themselves. In fact. 
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experience reveals that the toughest part of problem formulation 
is the prescription of proper boundary conditions. 

Surfaces (1), (4), (5), (9) have convective heat transfer 

boundary condition, with the overall heat transfer coefficient and 
ambient temperature given as input data. 

The weighted residual form of the convective boundary 
condition is - 


f dB = 


/( ^ 

§ N,h T - T ^ dB 

1 amb^ 


= - # N.h^T dB + # N.h'l dB (3.15) 

1 1 amb 

4 

Expanding the temperature as T = £ N .T . , the matrix 

contributions of convective heat loss boundary conditions are 
given by : 






dB 


= § N.h^N.dB tT.] - < 

^ J J J L 


N.h^T ^dB 

' I amb 


(3.16) 


On the surface(lO), heat flux is zero. Hence, 



(3.17) 


The rest of the surfaces are having ambient temperature 
boundary condition, and these surfaces are handled by supplying 
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boundary temperature value in the input data. 

In case of transient diffusion solution, only surfaces (2), 
(3), (7), (8) are having prescribed (given value) boundary 

conditions and other surfaces have zero flux condition. 

3.4 ELEMENT ASSEMBLY 


In the previous section, the elemental contributions to the 
left hand side coefficient matrix and right hand side vector have 
been discussed in detail for the four noded isoparametric 
quadrilateral elements. These elemental contributions are 
assembled into a global matrix equation by adding the entries 
corresponding to each nodal variable from neighbouring elements 
appropriately. 

The boundary integral contributions on both the left and 
right hand sides of the matrix equation are also incorporated as 
discussed in the previous section. This assembly procedure 
results in a matrix equation of the form : 

[A3{X} = {Bl 


where, 

[A3= coefficient matrix 

{X}= solution vector of temperature / concentration 
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obtain the following quantities ; 

(i) Isotherm plots 

(ii> Graphs for the variation of interfacial temperature 

(iii) Iso-concentration plots 

(iv) Iso- yield stress plots 

(v) Approximate wear profiles 

ISOTHERMS : Temperature has been predicted for all the nodes in 
the domain. - Using linear interpolation, it is possible to 
estimate temperature at any point across the element boundary. 
The coordinates of points with same temperature are calculated 
using this procedure and joined using smooth curves. At a glance, 
isotherms give a clear idea about the temperature distribution in 
the domain. 

ISOCONCENTRATION CURVES *. Once the temperature field is available 
in the vicinity of the cutting tool tip, the nodal values of 
carbon concentration can be calculated easily using equation 
(3.9). Iso-concentration curves are then plotted to visually 
depict the extent of diffusion transport, at <^fferent times. 
ISO-YIELD STRESS CURVES : Based on the nodal temperature and 
concentration values, the local yield stress has been calculated 
as discussd in chapter 2. Then contours of constant yield shear 
stress in the tool material have been plotted. Comparing this 
with the actual stress in the tool material, the wear profile can 


be obtained. 
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3.7 PROGRAM FOR FINITE ELEMENT ANALYSIS 

The programming work is divided into two parts : 

(i) Program for grid generation 
and, (ii) program for finite element analysis. 

In the grid generation program, coordinates of all the points 
within the solution domain and the elemental connectivity are 
calculated. These are used in the finite element analysis. 

Fig. (3.3) gives the flow chart of the analysis program. 

The main program calls the subroutines DINPUT and ITEEAT 
which in txirn call the other subroutines. 

Subroutine DINPUT reads and returns all the input data such 
as cutting conditions, initial conditions, boundary conditions, 
nodal coordinates, elemental connectivity data and prescribed 
parameters like ambient temperature and material properties at 
ambient temperature. The whole input data is in dimensional form. 
Some of the inputs are checked for correction by two routines 
named DIAGNl and DIAGN2. 

Subroutine ITERAT calls the FRONTS subroutine for solving 
assembled matrix equation. This subroutine begins the problem 
with initial values of all the nodal variables. This routine also 
calls WRITER subroutine to write the nodal values of the variables 
if the solution has converged. 
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Fig. 3.3 : Flow Chart of the FEU Program 
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FRONTS also calls the MATRIX subroutine which calculates the 
element fluid matrix 'fiumx'. This fluid matrix consists of the 
elemental contributions of the governing equations to form the 
coefficient matrix. Similarly elemental right hand side vector is 
also calculated here. Subroutine DRIVES is called within MATRIX, 
which in turn calls the subroutines SHAPE4 and DJACOB to calculate 
shape functions and their derivatives. MATRIX subroutine also 
cans subroutines VINTER, VCUT, HEATGEN and PROPERTIES for 
calculating quantities such as the inter facial velocity, cutting 
velocity, heat generation and material properties respectively for 
the temperature solution. Similarly, during the solution of the 
species diffusion problem, the MATRIX subroutine caUs the 
subroutines of CALDIFF, VINTER, VCUT and PROPERTIES for 
calculating the diffusion coefficient, inter-facial velocity, 
cutting velocity and material properties. 

The boundary conditions are implemented in BOUCON which is 
called in MATRIX before it returns to FRONTS. This is called for 
only boundary elements. 

3.8 CLOSURE 

Although the formulation and the programme code developed in 
the present study are general, they could not be applied to a 
variety of realistic situations involving the wear of HSS, WC aind 
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coated tools. The main reason for this was the lack of data on 
material properties, especially on the diffusion coefficients and 
and regarding the effect of alloying component concentration on 
the strength of the tool material. Also data on the machining 
conditions at high cutting speeds were scarce. Therefore, only a 
few model studies could be conducted and even here, diffusion wear 
was important only for some cases. With more complete data, the 
analysis can be applied to several practical situations to study 


tool wear. 


CHAPTER 4 


RESULTS AND DISCUSSIONS 


Using the Finite Element solution procedure discussed in the 
previous chapter, temperature and concentration fields have been 
predicted for a variety of typical machining conditions for HSS 
and carbide tools. These have been further processed to obtain 
the variations in the local yield stress value, finally leading to 
an estimate of the crater wear length in the case of the HSS 
tools. For carbide tools, the post-processing was not possible on 
account of data scarcity on variation of tool hardness. The 
results obtained from the analysis are discussed below. 


4.L TEMPERATURE RESULTS 

For HSS tools, resTilts of earlier researchers are available 
mainly in the low cutting speed conditions. Although the 
diffusion effect as well ais the tool wear will be negligibly small 
for such conditions, computations have been made for low speed 
situations also with the intention of comparing the temperature 
predictions. 

In Figs. {4.1-4.4), the steady-state isothermal contours 
within the toed, and the chip have been presented. The predicteed 
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TOOL MATERIAL = HSS 
a = 41'' 

V = 24.76 m/niin 



1. T = 826.285 K 

2. T = 786.285 K 

3. T = 746.285 K 
10. T = 466.285 K 
at = 40 K 

( IN TOOL ) 



1. T = 826.285 K 

2. T = 506.285 K 

3. T = 346.285 K 
AT = 80 K 

( IN CHIP ) 


Fig. 4.1 : Temperature Contours 
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TOOL MATERIAL = TUNGSTEN CARBIDE 
a = 20** 

V = 130.73 m/min 



1. T = 1178.28 K 

2. T = 1058.28 K 

3. T = 938.28 K 
AT = 120 K 

( IN TOOL ) 



1. T = 1178.28 K 

4. T = 818.28 K 
7. T = 458.28 K 
AT = 120 K 
( IN CHIP ) 


Fig. 4.2 : Temperature Contours 
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TOOL MATERIAL = HSS 
a = 25'' 

V = 140 m/min 



1. T = 1600 K 

2. T = 1460 K 

3. T = 1320 K 

4. T = 1180 K 

5. T = 1040 K 


( TOOL ) 


Fig. 4.3 : Temperature Contours 
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TOOL MATERIAL = HSS 
a = 20^" 

V = 154.6 m/min 



1. T 

2. T 

3. T 
AT = 


= 1037.18 K 
= 897.18 K 
= 757.18 K 
140 K 


Fig. 4.4 : Temperature Contours 
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contours are in reasonable agreement with the restilts of Murarka 
et al. C19793 and Tay et al. El 9743. The maximum temperature 
isotherm occurs slightly away from the cutting edge on the 
chip-tool interface. This is obviously the result of metal 
getting heated up through primary and secondary deformation zones 
sequentially. In the chip also, high temperature zone occurs 
adjacent to the interface as expected. It is clear that the 
maximum temperature increases with cutting speed, due to the 
larger rate of heat production {caused by the plastic 
deformation). 

The variation of interfacial temperature with distance from 
the cutting edge are shown in Figs. (4.5) and (4.6), for two 
different cutting conditions. The predictions of the present 

analysis have been compared with those of Murarka et al. [19793. 
Although there is reasonable agreement in the predicted trends, 
some differences are also observed. In the present analysis, the 
variations in material thermal properties such as the thermal 
conductivity and specific heat have been fully taken into account. 
Also a variable shear stress has been assumed to exist on the 
chip-tool interface as known from independent stress measurements. 
Murarka and co-workers assumed constant properties and uniform 
interfacial shear stress. Their predictions were justified on the 
basis of comparing the theoretical resxilts for the average tool 
tip temperature with actual measurements by tool thermocouple 
method. The average tool tip temperature predicted from the 





TEMPERRTURE (K) 


Chapter 4 •• Results and Discussions 


70 


TOOL MATERIAL = HSS 

MURARKA et al. 

PRESENT MODEL 


T { MURARKA ) = 718.76 K 

av 

a = 41® T ( PRESENT MODEL ) = 730.45 K 

V = 24.76 m/min 



Fig. 4.6 ; Variation of Interfacial Temperature 
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present analysis is quite close to that obtained by Murarka et al. 
El 979]. 

4.2 CONCENTRATION PROFILES 

In Fig. (4.7), the concentration profiles at the various 
times are presented for the carbide tool. It is evident from 
these figures that diffusion of carbon within the tool gets 
depleted in course of time in the region close to the interface. 
The depth affected by carbon diffusion increases with time in the 
beginning, but eventually a steady-state is approached. The 
reason for this lies in the fact that the high temperature zone 
occurs close to the chip-tool interface. Diffusion being a 
process which is very sensitive to the temperatutre (exponential 
dependence), a region close to the interface is affected. The 
concentration profiles obtained for the HSS tool (Fig. (4.8)) are 
similar to those for the carbide tool. 

A point to be noted regarding the results of Figs. {4.7-4.8) 
is that these correspond to fairly high cutting speeds. At low 
speeds, the temperature level is also low and diffusion is 
virtually non-existent. 

4.3 YIELD STRESS CONTOURS 


In Figs. (4.9 and 4.10), the constant yield stress contours 
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TOOL MATERIAL = TUNGSTEN CARBIDE 
a = 20® 

V = 130.76 m/min 



1. c = 0.557 % 

9. c = 0.621 % 

10. c = 0.629 % 

Ac = 0.008 % 
AFTER 5 SECONDS 



1. c = 0.295 % 

8. c = 0.645 % 

9. c = 0.695 % 

Ac = 0.05 % 

AFTER 15 SECONDS 



1. c = 0.174 % 

8. c = 0-734 % 

9. c = 0.814 % ■ 

Ac = 0.08 % 

AFTER 25 SECONDS 


Fig. 4.7 ; Contours of Carbon Concentration in Tool 

Al 
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TOOL ^TERIAL = HSS 
a = 22 

V = 150 m/min 
TIME = 10 min 



1. a = 501.714 N/mm^ 

I o 

2. a 541.714 N/mm^ 

3. 0^^= 581.714 N/mm^ 

9. 821.714 N/mra^ 

Acr = 40 N/mm^ 

1 tj' 


Fig. 4.9 : Yield Stress Variation in Tool 
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TOOL MATERIAL = HSS 
ce « 20® 

V = 154.6 m/nsin 
TIME » 10 idn 



1. a = 501.69 N/mm^ 

YS 

2. o = 601.69 N/mm 

2 

3. cr„ = 701.69 N/mm 

6. a = 1001.69 N/mm' 
2 

A<7 = 100 N/mm 

YS 


Fig. 4.10 : Yield Stress Variation in Tool 
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Within the tool have been plotted. The material strength 
decreases considerably (almost by a factor of 5-6) near the 
chip-tocil interface due to both high temperature conditions and 
the depletion of the alloying element (carbon). The contours have 
shapes approximately similar to that of the crater. If a detailed 
stress analysis is done for the tool side, the yield stress 
contours will be useful for applying local yield criterion to 
predict the depth of crater wear. It must, however be cautioned 
that wear being a dynamic phenomenon, a dynamic stress analysis is 
necessairy. 

4.4 CRATER LENGTH PREDICHONS 

In Figs. (4.11-4.15), the yield stress values along the 
interface are compared against the assumed actual stress 
variation. Since it is reasonable to expect that wear should 
occur if yield stress is less than the actual stress, the zone 
where this ineguality holds should be prone to wear. Thus, the 
approximate length of the crater can be predicted for a given 
maohlnina time by oompartaa the local yield stress and actual 

stress value. 

in the present warb, due to time UudtaUons, it was not 
possible to perform a complete stress analysis o£ the tool 
portion, coupltog the stress predictions with transient thaimal 
and aiffudon analysis vdthln the tool, the extent of diffusion 
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TOOL MATERIAL - HSS 
TIME - 10 min 
“ - - YIELD STRESS 

ACTUAL STRESS 


a = 20'' 

V = 120 m/min 
X = 1634.1 K 

max 

fh = 3587.7 N 
fv = 1190.7 N 



Fig, 4.11 : Comparison of Interfacial Stresses and Estimation 

OF Crater Length 
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TOOL MATERIAL - HSS 
TIME - 10 min 
~ YIELD STRESS 

ACTUAL STRESS 


a = 25° 

V = 140 m/min 
T = 1622.81 K 

max 

fh = 3610.2 N 
fv = 1210.8 N 



Fig. 4.12 : Comparison of Interfacial Stresses and Estimation 

OF Crater Length 
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TOOL HATERIAL « HSS 
TIME » 10 min 
" ■ YIELD STRESS 

■* - ACTUAL STRESS 


« = 22 ^ 

V = 150 m/rain 
T = 1628.96 K 

max 

fh = 3604.6 N 
fv = 1200.3 N 



DISTRNCE JTO1 CUTTING EDGECmm) 


Fig. 4.13 ; Comparison of Interfacial Stresses and Estimation 

OF Crater Length 
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TOOL MATERIAL - HSS 
TIME « 10 min 
~ - ~ YIELD STRESS 

-- - ACTUAL STRESS 


a = 20® 

V “ 154.6 m/min 
T = 1037.18 K 

max 

fh = 3500.0 N 
fv = 1560.0 N 



Fig. 4-14 : CoHPAiis©ii of Ihterfacial Stresses and Estimation 

m Crater Ls^gth 
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TOOL MATERIAL - HSS 
TIME - 10 min 
- ■ - YIELD STRESS 

ACTUAL STRESS 


cc = 20*" 

V = 160 m/min 
1^5L67 K 

fh = 3520 N 
fv = 1260 N 



Fig. 4.15 : Compariscw of Ihtehfacial Stresses and Estimation 

OF Crater Length 
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wear at any instant can be judged. The results of the present 
study indicate that diffusion wear is important only at elevated 
temperature conditions (above 900*^ C). This observation is in 
agreement with the trends known from practical measurements. 


4.5 CONCLUSION 

A Theoretical model for estimating the diffusion wear has 
been proposed, based on the temperature and concentration fields 
within the tool region. The model proposed here considers carbon 
diffusion only. However, diffusion of other species can also be 
incorporated in a similar way and the effect of such diffusion 
phenomena on tool strength can be analysed. Finally, based on a 
detailed stress analysis, it is possible to obtain wear profiles 


at different times. 



CHAPTER 5 


SUGGESTIONS FOR FUTURE WORK 


In the present work, only the heat exchange and diffusion 
equations have been solved for the sake of simplicity. These can 
be coupled with a visco-plastic model for the chip region and an 
elasto-plastic stress analysis of the tool region. 

Diffusion of other species such as tungsten, cobalt, 
chromium, vanadium, molebdenum etc. can be considered. The 
effects of such difussing species on the tool strength can be 
studied. 

Considering a blunt tool, flank wear can be studied in a 
similar manner. 

Different cutting tool materials {various HSS grades, carbide 
tools, coated carbide tools etc.) can be studied. 
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